Zero-momentum modes and chiral limit in compact lattice QED * 
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The influence of zero-momentum gauge modes on physical observables is investigated for compact lattice QED 
with dynamical and quenched Wilson fermions. Within the Coulomb phase, zero-momentum modes are shown 
to hide the critical behavior of gauge invariant fermion observables near the chiral limit. Methods for eliminating 
zero-momentum modes are discussed. 



Lattice gauge theories allow to compute physi- 
cal observables without gauge fixing. At the same 
time, more detailed information about nonper- 
turbative properties of quantum fields can be ex- 
tracted from gauge-dependent objects. But the 
straightforward application of iterative gauge fix- 
ing procedures leads to the appearance of gauge 
(Gribov) copies. For compact lattice QED in the 
Coulomb phase, such spurious gauge configura- 
tions disturb the correct perturbative behavior of 
both the photon and fermion correlators 
Numerical |^-|^ and analytical Q investigations 
of this problem have shown the main excitations 
responsible for the occurence of such copies to be 
double Dirac sheets and zero-momentum modes. 
The former can be removed by usual, or periodic, 
gauge rotations. However, the removal of the 
latter requires constant, or nonperiodic, gauge 
transformations, which - for finite lattice size - 
violate the invariance of the fermion matrix. 

In this talk we present results of a study of 
the influence of zero-momentum modes in com- 
pact lattice QED within the physical (Coulomb) 
phase. 

The action S is given by 

S = SG[U] + SF[U,i>,i;], (1) 
Sg = /3 ^ (l~ReC/,^C/,+^,,C/,V,^C/:,), 
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Sf = EE^'-^-v^y' 

f=l x,y 

where A4 is the Wilson fermion matrix 

Mxy = ^D^xy - '^Ej^-'-I^ ^ lti)Uxti5y,x+[i (2) 

+ (1d +7^)^^^?/^^-*^} • 

Ux^l, = cxp(i6'a;^), 9x^l, £ (-7r,7r] denote the 
link variables. (3 and k are the inverse square 
bare coupling and the hopping parameter, respec- 
tively. The gauge as well as the fermion field 
obey periodic boundary conditions (b.c.) except 
for the X/^ (time) direction, for which the fermion 
field will be taken either periodic or antiperiodic. 

In what follows we consider both the quenched 
case simulated with a heat bath method and the 
dynamical case for A^f = 2 flavors studied with 
the hybrid Monte-Carlo method (HMC). 

When studying gauge-dependent objects like 
gauge link correlators or the fermion propagator 
the Lorentz (or Landau) gauge fixing procedure 
is usually applied by iteratively maximizing the 
functional 

^["] ^ E ^^^^ ' (3) 

with U% = n.,U^^W^+^ and e U(l). V - 
iVg • iV4 denotes the lattice volume. The gauge 
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Figure 1. Pion norm as function of k for full (and 
quenched) compact QED with Wilson action at 
(3 — \.l for various lattice sizes (data taken from 
§)• 



functional in most cases has many local extrema. 
Consequently the iterative maximization proce- 
dure can provide different gauge copies. But the 
'best' or physically correct gauge copies are ex- 
pected to be related to the global maxima. 

In our recent paper 0| gauge variant fermion 
correlators in the quenched approximation have 
been calculated within the standard Lorentz 
gauge. It has been shown that these correlators 
behave strongly decaying, as if they were mas- 
sive in the limit k — > Kc(/3)- This effect has been 
shown caused by constant or zero-momentum 
modes of the gauge fields 



(4) 



Within the Lorentz gauge the zero-momentum 
modes do not allow to obtain the best gauge copy, 
i.e. the maximization of the functional F[^] de- 



fined in eq. (||) does not lead to its absolute max- 
imum. Zero-momentum modes cannot be elimi- 
nated by usual (periodic) gauge transformations. 
The Lorentz gauge can only drive them into the 
interval [— tt/TV^, -I-tt/A^^] which, however, ap- 
pears to be sufficient to distort gauge-dependent 
observables. Increasing the lattice volume does 
not change the situation significantly. 

In Q an iterative Lorentz gauge fixing proce- 
dure combined with zero-momentum mode sub- 
traction was proposed (ZML gauge). The zero- 
momentum mode suppression 



z^,[U] = 



(5) 



is achieved by appropriate non-periodic gauge 
transformation steps 



eU(i) 



(6) 



following each Lorentz gauge step. 

Within the ZML gauge the Lorentz functional 
(^ reaches its global maximum with high accu- 
racy |^,||. On the other hand, ZML gauge fix- 
ing leads to the correct perturbative behavior 
of gauge-dependent objects In particular, 

the results are compatible with vanishing fermion 
masses near the chiral limit. 

In this talk, instead of gauge variant objects 
we consider the gauge invariant - with respect to 
periodic gauge transformations - 'pion norm' 



(n) = ( ^Tr 



where (■ ■ •) means the functional average with re- 
spect to the compact U{1) gauge field variables. 
The Ai are the eigenvalues of 75A^. One expects 
(n) to be a good indicator of the chiral limit at 
K as some of the Aj are expected to tend 

to zero. It is worth noting that for periodic and 
time-antiperiodic b.c. for the fermion fields the 
averages (11) strictly coincide. 

However, the numerical study of fermionic ob- 
servables like (n) near the chiral limit does not re- 
veal the critical properties as expected from low- 
est order and finite lattice size perturbation the- 
ory. This can be seen from the K-dependence of 
the pion norm numerically computed at low /?- 
values within the Coulomb phase ||] (see Fig. |^). 
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Figure 2. Pion norm in the free fermion case: 
without zero-momentum modes and in the zero- 
momentum mode background, lattice size 6^, pe- 
riodic (p.b.c), as well as time-antiperiodic (t- 
a.b.c.) boundary conditions. 



Figure 3. Pion norm as function of k, redefined 
with ZML gauge ((H)') for full and quenched 
QED, as well as without any gauge fixing ((11)) 
for fuU QED; all data for j3 = 2.0, lattice size 4'', 
periodic b.c.'s. 



Its behavior is very smooth and no sign of any 
critical behavior is observed. The volume depen- 
dence of (n) is rather weak, and there is no sig- 
nificant difference between the quenched and the 
dynamical case. 

It is interesting to compare these results for (11) 
with the free fermion case given by 

(n)o = ^HWH^^^'-W"^ (7) 

p \ i-i ^ 

-^(--i:-^)y\ 

where the p^, /i = 1 , • • • , 4 are integers except for 
time-antiperiodic b.c. causing p4 to take half- 



integer values. In Fig. || one can see the k- 
dependence of (11) o calculated on a symmetric 
lattice (A'4 = A^s = 6) for periodic and time- 
antiperiodic b.c. For periodic b.c. (n)o obvi- 
ously gets singular at k = 1/8, whereas for time- 
antiperiodic b.c. the K-dependence of (n)o be- 
comes smooth for symmetric lattices. However, 
in the latter b.c. case (n)o develops a peak 
for strongly elongated lattices {N^ — » 00 with 
A^,=fixed). 

We are going to demonstrate that this drastic 
difference between (H) (Fig. |l|) and (n)o (Fig. |) 
is due to constant or zero-momentum modes of 
the gauge fields Zfi\U] as defined in eq.(^. 

For the free case {j3 00) this can be easily 
demonstrated by integrating (11) over constant 
modes. Both for periodic and time-antiperiodic 
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b.c. it yields 

+ (i-2k^cos<?!>^)U . 

The latter expression is completely smooth in k 
and agrees with the former time-antiperiodic, free 
result for symmetric lattices (see Fig. ||). 

Now let us consider the interacting case, i.e. 
for finite /3-values. We redefine the pion norm 
such that zero-momentum modes become elimi- 
nated. Each gauge field configuration generated 
in the simulation is transformed subsequently by 
the ZML gauge procedure including non-periodic 
gauge transformations in order to suppress zero- 
momentum modes as described above. We aver- 
age n with respect to the ensemble of transformed 
fields. Note that the fermionic part of the action 
is not invariant under constant gauge transforma- 
tions @. Therefore, we get a new average (11)' 
which differs from (11). 

In the following we choose periodic boundary 
conditions, because we expect from the free case 
that they lead to a more pronounced chiral be- 
haviour than the time-antiperiodic ones. 

In Figure ^ we show the dependence of the pion 
norm (11) and (11)' on k. One can see that for 
dynamical fermions (full circles) as well as for 
quenched fermions (boxes) the redefined observ- 
able (n)' has a sharp singularity near the point 
Kc = 0.1307(1) for P = 2.0 §. In contrast, the 
standard definition of the pion norm (11) demon- 
strates a completely smooth behavior (open cir- 
cles). We checked these results for (11)' also on 
larger lattices. For n approaching Kc the same 
critical behaviour is observed, whereas very close 
to and slightly above Kc the infiuence of an in- 
creasing number of very small fermionic eigen- 
modes leads to stronger fluctuations ('exceptional 
configurations'). The dynamical and quenched 
results resemble each other. This can be inter- 
preted such that the zero-momentum modes - al- 
though removed from the observable (11)' - con- 
tinue to dominate the fermionic determinant. 

What about other (maybe simpler) methods to 



get rid off the constant modes @) of the gauge 
fields? We have also considered the Polyakov 
line gauge invented in ^ which transforms the 
spatially averaged Polyakov line values into real 
numbers. We convinced ourselves that this non- 
periodic gauge - without the necessity to employ 
the Lorentz gauge - leads to the same singular 
chiral behaviour of the pion norm as the ZML 
gauge. 

Our main conclusion is that zero-momentum 
modes play an important role near the chiral limit 
of compact QED in the Coulomb phase, at least 
for finite lattices. They smooth out the criti- 
cal chiral behaviour expected from lowest order 
perturbation theory for observables like the pion 
norm and the chiral condensate. In any case it 
is necessary to take the zero-momentum modes 
properly into account. 
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